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STEREO ELECTRO-OPTICAL TRACKER STUDY 
FOR THE MEASUREMENT OF 

MODEL DEFORMATIONS AT THE NATIONAL TRANSONIC FACILITY 

RICHARD J. HERTEL AND KENT A. HOILMAN 

( 

SUMMARY 


Model deformations must be measured if aerodynamic 
parameters are to be accurately determined. A report^ recently 
completed investigated the feasibility of using a stereo tracker 
system to measure deformations in the National Transonic Facility. 
The present study investigates the influence on the stereo tracker 
system by the unique conditions found in the NTF. None of the 
conditions examined, (i.e., shock waves, boundary layers, 
pressure, temperature, vibration, and optics) will prevent 
meaningful deformation measurement provided care is taken in 
system design and the data is post-processed. 

The use of a CID type system was examined and found to 
be inadequate for deformation measurement. Even with the expected 
technology advances, the proposed image dissector camera will 
remain superior in the near future. 

A computer program was developed to trace virtually 

any ray. through the co_mple_x _ tunnel conditions . The program will 

handle shock waves, boundary layers, pressures, temperatures, 
nonhomogeneous index of refraction, and optics of any conic sec- 
tion or polynomial of sixth order or less. This program 
approaches the goal of a universal ray trace program capable of 
handling any tunnel condition or ray geometry. 
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INTRODUCTION 


The primary purpose of this study was to investigate 
the effects on measurement accuracy of vibration, camera and win- 
dow non-linearities and potential disturbances in the system opti- 
cal path. Therefore, a detailed understanding of the system is 
not required but a familiarity with the general data gathering 
process will be helpful. The Stereo Electro-Optical Tracking 
System (SETS) measures the position of light emitting targets on a 
wind tunnel model by tracking the targets with two or more image 
dissector cameras located in the floor of the wind tunnel. These 
cameras or trackers are capable, under computer control, of measuring 
the brightness of a selected position of the scene. By measuring 
scene brightness of several areas the location of the targets can 
be determined. The exact target location is found by an iterative 
process which converges on the center of the target image and does 
no- rely on a single measurement or on a well defined target 
image. The location or positions of targets as measured by two or 
more trackers are combined in a stereo algorithm to determine true 
target position. Each tracker is capable of locating 1700 targets 
per second. Table 1 lists some of the SETS characteristics. 

Figure 1 shows a typical tracker configuration and its mounting 
arrangment. The inside of the tracker enclosure will be main- 
tained at room temperature and atmospheric pressure. 

Figure 2, 3, and 4 show the location in the tunnel 
for the six possible trackers. All of the trackers will be 
mounted below the model. All targets on the model are viewed by 
at least two trackers. Small luminous targets will be built into 
the model and controlled by the SETS electronics. Each target 
will be activated only when the trackers are viewing it. The 
entire system will be under computer control thus allowing the 
targets to be measured many times in any order. 
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Table 1 


SETS 


Sensed Volume 


Target Sources, Number 
, Type 
, Size 


Model Deformation at 
Target Location 

Sensor Position Deviation 
Detection 


Characteristics 

0.8m by 1.4m by 0.5m 
5 to 5 0 

Light-Emitting Diode 
0.5mm diameter 


+ 20mm maximum 


0.05mm all axes 


Model Vibration Limits 


+ 12.5mm at 20 Hz 
+ 5mm at 50 Hz 
+ 1.25mm at 200 Hz 


Target Sampling Interval 
(Single Target Location) 

Target Sampling Rate 

Sensor Tube Type 

Sensor Environment Control 

Tracking Electronics 
Calibration, geometric 
Pre-Test Reference 
Data Processor Type 
Data Storage 
Control/Display 
Hard Copy 


590 microseconds, typical 

1695 targets per second, typical 

ITT F4012, S20 photocathode 

70°F Dry Air at 130 psi (NASA 
Supplied) Electrical Heating 

Digital, in Tracker head 

Area map, prior to test 

Targets on Stable Portion of model 

MODCOMP 

Magnetic Tape 

CRT Display/Keyboard 

Line Printer/Plotter 
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FIGURE 1 TUNNEL AXES AND TRACKER POSITION 
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FIGURE 2 TRACKER AND FIELD LOCATIONS 
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FIGURE 3 TRACKER COVERAGE OF MODEL SPACE IN XZ PLANE 



QUARTZ WINDOW 



FIGURE 4- MULTI-TARGET TRACKER ASSEMBLY 
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ERROR SOURCES 


Successful operation of the SETS requires that the 
trackers accurately measure the position of each model target in 
the wind tunnel environment which includes vibration of targets 
and trackers, extremely low temperatures, high static pressures, 
and potentially significant disturbances in the optical path of 
the cameras which result naturally from operation of the wind tun- 
nel. 

Eight separate areas of concern have been identified 
and are examined in this section. Some areas prove to be no 
problem. Some simply require the correct system design to mini- 
mize them. Some require that the data be corrected in a post pro- 
cessing step. All of the areas are manageable and should allow 
successful operation. This is summarized in Table 2. 


Table 2. Error Summary 


Area Investigated 

Action Necessary I 

Error Source 

No 

Problem 

System 

Design 

Correction 

Required 

Window Effects 

X 

X (tilt) 


Tunnel Pressure 
& Temperature 



X 

Laminar Boundary 
Layer Effects 

X 



- - -Turbid Boundary — - . 

Layers 

. X _ _ 

- _ X. (Target . 
Size ) 


Shock Wave 
Effects 



X* 

Vibration 



X* 

Targets 


X 



*System Performance Limiting Effect. 
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Tunnel Window Effects 


Window effects upon the optical measurement are due to 
the abrupt difference in index of refraction between the window 
and tunnel atmosphere. This difference coupled with the finite 
thickness of the window, the window curvature due to temperature 
and pressure, and the off-normal incidence of the camera line of 
sight cause a deviation in the optical path between the camera and 
the model. 


Window Focus Effects 

The curvature of the window by temperature or pressure 
introduces a small amount of additional optical power, Ap , to the 
system. 

Ap = (n - l)/ri + (1 - n)/r 2 - t(n - 1) (1 - n)/(n r^ r 2 ) 
where n = index of refraction of the window, quartz 1.4581 

t = thickness of the window 
rj^ = inner radius of curvature 
r 2 = ti + 1 = outer radius of curvature 
When the window is concave on the model side, r^^ and r 2 are posi- 
tive valued. 


According to the data provided by NASA, the radius of 
_ curvature of a . 1-inch thick window having a temperature gradient 
of +85°F across its faces is 37950.66414 inches. (Center sag is 
110 microinches.) 


Ap = (1.4581 - 1. )/37950. 66414 
+ (1. - 1.458D/37951. 66414 

- (1)(0.4581)(-0.4581)/(1.4581)(37950.)(37951.) 

= 4.18 E-10 inches"! (1.65 E-11 mm"!) 

A p is the reciprocal of the focal length of the window 

"lens". 

1/Ap = 2.39 E+9 inches (6.08 E+lOmm) 
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This is a very weak lens. The camera optics are 
approximately 16mm focal length or Ap = 0.06mm"l. Temperature 
deformation of the window does not appear to be a significant fac- 
tor affecting focus. 

Pressure changes also cause the window to sag. The 
sag is 33 microinches, typical, to 54 microinches maximum 
according to NASA data. These sags are one-half or one-third 
those introduced by the temperature gradient. Again no signifi- 
cant focus effects are expected. 

Window Displacement Effects 

Figure 5 shows a ray incident on the window at an 
angle 0 . This angle causes the ray path to be displaced from a 
straight line by an amount d. 

d = t sin 0 / r 

V (1 - sin 0) 


■^1 

, 2 2 . 2 
^2 ~ 


) 


If the window tilts, 0 changes by a small amount d O. 
The sensitivity of d to changes in 0 is : 


^ d0 = -i- 

90 - 30 


t sin 0 


2 h 

a - sin 0) 


, 2 2, . 

(n^ - n^ sin 0) 


)) 


d0 


= 19.2 mm/radian at 0 

^1 

t 


= 35° 

= 1.00298 
= 25 mm 


The objective of the measurement process is to locate 
the targets within 150 micrometers uncertainty. Therefore, the 
maximum tilt during any series of target measurements is: 
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FIGURE 5j Refraction at Test Section Window 
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A0 
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150 um „ 

-— radian = 7.8 mrad 

19.2mm 


There is no reasonable estimate of the likely tilts 
available at this time. 

The window must also have parallel faces. If the win- 
dow has a wedge cross section, displacements of the window in the 
plane of the tunnel wall will cause changes in d. 

||- dt = 0.27 dt 0 = 35° 

n^ = 1.00298 
U2 = 1.4581 

The maximum thickness change must be less than: 

Atm < < — = 550 ym 

0.27 


over the maximum likely lateral displacement range. If 10 mm is 
used as the largest displacement range, then the maximum wedge 
angle is an easily-achieved value: 


Tunnel Pressure and Temperature 



During set up and initial target acquisition, the tun- 
nel will be at normal atmospheric pressure and room temperature.- 
During actual data runs, the tunnel may be cooled and pressurized. 
Both cooling and pressurizing the gas in the tunnel increases its 
density and thereby increases its optical index of refraction. 

This increase in the medium surrounding the model makes the model 
appear larger and distorted. 
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The index of refraction of nitrogen varies with the 

temperature and pressure by the following formula; 

(n - 1) 

. ^ o P 

(1 + 0.00366 t) 760 

where Uq = 1.000298 for 0°C and 760 mm Hg. 
thus n(-195°C, 9 atm) = 1.00937 

n(+95°C, 1 atm) = 1.00022 

If the tracker is maintained at a constant pressure 
and temperature while the tunnel is pressurized and cooled, the 
maximum image change possible is produced. The interaction at the 
tunnel window can be reduced for calculation purposes to a single 
interface where Snell's law (nsin© = n'sin© ') applies. Using 
the tunnel geometry and a FOV of 40“ for the tracker, the effects 
can be calculated. Figure 6 shows the situation for a tracker 
located at position 4. The 40 degree FOV allows a model 1600 mm 
long to be viewed. 

Table 3 gives the angles and dimensions as viewed from 
the tracker for the two extreme tunnel conditions. The model 
appears to grow by approximately 4 percent. This scale change is 
not linear. The scale change of 4 percent can be followed by the 
tracker as discussed in the next paragraph; however, the non- 
linearity introduced must be removed by the data analysis 
- - routines . 


If the tunnel conditions change at such a rate that 
the target with the maximum optical induced motion is just 
tracked, then that rate is the upper limit on change allowable for 
the gas temperature and pressure. If the capture range of the 
tracker is a 200 x 200 ym square and the target is viewed 30 times 
per second, the maximum motion is then 5600 y m/s on the tube face 
or (using a 250 x scale factor) approximately 350 mm/s on the 
model. Thus, the temperature and pressure must change from room 


13 



FIGURE 6 TYPICAL OPTICAL GEOMETRY 
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0 Before 
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£ After 
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Before Conditions: +95C, 1 atm 

After Conditions: -195C, 9 atm 

See Figure 7 for angle and length definitions 


A£ Change 


4.4 

8.8 

16.5 
23.2 
33.7 

73.6 
















/ , 


n 


ambient to LN 2 and 9 atmospheres in less than 1/4 second to cause 
loss of track. This will not be an operational limitation, as 
cooling requires hours. 

In order to calculate more accurately the influence of 
pressure and temperature changes, the simple optical model of 
Figure 6 was modified to that in Figure The tracker is assumed 

to be on the tunnel centerline (position no. 4 of Figure .2) and 
tilted at an angle of 50® in the X-Z plane. The 50° angle provi- 
des coverage of the entire model with a reasonably well-centered 
image. The tracker optics are positioned 20 mm outside the tunnel 
window and consist of an ideal 16 mm lens. 

The model of Figure 7 was programmed in Fortran IV. 
Program ITO traces rays from the image in the tracker to the 
object in the tunnel. Program OTI traces rays from the object in 
the tunnel to the image in the tracker. This ray cannot be traced 
directly and is found by an iteration technique. Tables ^4 and S 
are sample printouts from the OTI program. Tables .6' and 1 ' are 
sample printouts from the ITO program. The scale is the size of 
the output variable divided by the size of the input variable. 

Thus for the ITO program, the scale is the size of the image 
divided by the size of the model for the specific image and object 
conditions used. The scale varies from 0.013 (the image is 1.3 
percent of the size of the model at the tail to 0.004 at the nose). 
This scale is in the x axis of the tunnel only. For the OTI 
program, the scale is the object size divided by the image size. 

If a particular ray is traced in both directions, the object and 
image points should match and the scales will be reciprocals. 

When the tunnel is pressurized to 9 atmospheres and 
cooled to LN 2 temperatures, the increased density of the gas makes 
the model appear approximately 4 percent larger. The size 
increase is not uniform or linear over the model. The tail grows 
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Camera Optical Axis — 50*^ 

TAN a = Z/F 
TAN Y = X/H^ 
x \2 Sin e = Sin y 
a + e = 40° 

Tan B = R/S 

Figure ,7. Simple 


= 1200 mm 
S = 20 mm 

H 2 = 1200 mm 
F = 16 mm 

n^ = 1.00937 (9 atm, LN 2 temp) 
n 2 = 1.00022 (1 atm, Room temp) 

)tical Model 
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Table 4 

OBJECT TO IMAGE RAY TRACE 
Room Tem perature , One Atmosphere 
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Table 5 

OBJECT TO IMAGE RAY TRACE 
LN 2 Temperature, 9 Atmospheres 
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Table 6 
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about 1 percent while the nose grows over 9 percent. This change 
becomes important when the model's absolute position is calcu- 
lated. The apparent shifts caused by the density change could 
possibly be interpreted as model motion caused by air loadings or 
structural changes. The local scale factor between object and 
image also increases with gas density. The scale increases vary 
from 1 percent at the tail to 9 percent at the nose. The scale is 
important when vibrations are being considered or when measurement 
precision is determined. The absolute position of the target is 
not necessary if a vibration measurement is made. Only changes in 
its position relative to some reference is needed. If the model 
grows, the reference position simply shifts. If the scale changes, 
the vibration measurement changes proportionately. 

The distortions introduced by pressure and temperature 
must be removed in the analysis of the data. The correction can 
be calculated provided the tunnel conditions are known. 


The optical model of Figure 7 can be described by the 

sz 

following equation: Q S tan C - — 

^'7 ' ^2 “”” 


\ 




1 + — tan C 


where , Q 


n^ F sin C 


Z cos C 


/2 2 
n„ /Z + F 


/2 2 
n^ /Z + F 


C = angle of camera 
Z = image location 
D = target location 
F = focal length 
Hi = tunnel height 
H 2 = tunnel length 
ni = index in tunnel 
n 2 = index outside tunnel 

The scale number is then the derivative of this 
equation and is: 


/ 

s 
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dZ 


dZ 


S/F 







1 



C 


s 

-= S tan C) tan C 

r 

( 1+1 tan C)^ F 

r 


where ; 


dz 


- cos C 

/2 2 
y'Z + F 


Z (F sin C - Z cos C) 

n^ (Z^ + F^) v4^ + F^ 


Laminar Boundary Layers 

Boundary layers will be present on the wind tunnel 
test section walls and on the surface of the model. These layers 
have the effect of changing the angle of rays passing through them 
and of displacing the rays. This can be demonstrated by approxi- 
mating a boundary layer as a series of uniform layers. 



Figure 8. Layer Model for Boundary Layers 
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For each interface between layers Snell's law applies: 
nj^ sin 0 1 sin 02 

Thus for a series of layers: 

n]_ sin 0]^ = n2 sin 02 = n3 sin 03 = n4 sin 04 = etc. 
For a series of parallel layers, the product of the index of 
refraction and the sine of the angle is constant. Thus, the con- 
ditions of the first and last layer are all that are needed to 
compute the angle. Figure 8 shows the typical ray path through a 
boundary layer. 



Figure 9 . Light Path through Boundary Layer 
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The ray emerges from the boundary layer at a different 
angle and the light source appears to be shifted from its actual 
location. These two effects, angle and shift, must be handled 
separately. 


When viewing a light source immersed in a boundary 
layer, the source appears shifted as shown in Figure 9. For a 
single tracker, the shift is along the surface, but for a pair of 
trackers and a stereo transform, the shift is also into the sur- 
face. Either component of the shift depends on the ray angle, the 
boundary layer thickness, and the index gradient in the boundary 

layer. The index can be computed using the Gladstone - Dale 
2 > ‘3 

constant K: ' ■ 

n - 1 = Kp; where p is density 


The density for boundary layers can be related to tem- 
perature by:^'^\ p^ T, 


— = — 
Pi - 


The temperature can be found by: 


.' 4 , 


^ = /I + e _/ 

1 ^ 


Stine, H: Application of Scattering Theory to the Measurement 

. of Turbulent Density Fluctuations by ah~ Optical Method. MACA 
Technical Note 3719. 

3 - . 

' Liepmann, H.W.: Deflection and .^ffusion of a Light Ray Passing 

Through a Boundary Layer. Douglas Aircraft Report No. SM-14397, 
1952. 

4' 

'' Dommash, D. : Airplane Aerodynamics. Pitman Publishing Corp., 
1967. 
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where e is the boundary layer recovery factor. Thus, for a worst- 
case condition there is a 25 percent change in the quantity n - 1. 



.where M = 1.2 (Tunnel maximum)^ 

e = 1.0 (Total Recovery) 

Y = 1.47 (Nitrogen at -181®C) 


by 


The thickness of the boundary layer can be calculated 
5.2 X 


t = 




where = 120 x 10^ (tunnel maximum) 

X = 1 meter, (typical model fuselage length) so that: 
5200 


t = 


10954 


= 0.47 mm 


A worst-case situation can now be calculated by 
assuming the maximum index variation takes place at the free 
stream edge of the boundary layer and the entire boundary layer is 
at the lower index value. 


The shift in apparent target position can now be 
calculated as shown in figure 10. The shift along the surface is; 


^ Dommasch, D. , Airplane Aerodynamics. Pitman Publishing Corp., 
1967. 

^ Nickes, O.W. ; Status and Operational Characteristics of the 
National Transonic Facility. AIAA 10th Aerodynamic Testing 
Conference paper, 1978. 

6 Hansen, A.G.: Fluid Mechanics. John Wiley and Sons, Inc., 1967 

^ Shames, I;H.: Mechanics of Fluids. McGraw-Hill, 1962. 
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FIGURE ID Target Shifts due to Boundary Layer 




h L>/ 1 1 


= t ( tan 0 I 
Sg = t ( tan (j) 


tan e 2 ) 
tan (j) 2 ) 


(Camera A) 
(Camera B) 


The shift seen by a stereo pair is: 


s = 

H 


X + 1 Y + 1 


(horizontal along surface) 


£ cos (j)^ cos 0^ 

s = — : — — — (vertical into surface) 

V sin((J)2 + 02> 


where 

a 

= t ( tan 0 jL + 

tan 


X 

= tan 02 /tan 

<j)2 


Y 

= tan 0]^/tan 

h 

n]^ sin 

01 

= n 2 sin 02 


n 2 ^ sin 

<i>i 

= n 2 sin 




20° < 02 < 

60 



V 

V 

0 

0 

<N 

60 


nj^ = 1.00700 (maximum index in boundary layer) 
n 2 = 1.00937 (maximum index in tunnel) 


I 


This produces a worst-case shift of 7^<ym along the 
model when viewed by one tracker. When viewed by two cameras, the 
worst-case condition is a shift into the surface of 4.5 ym and 
1.0 ym along the surface. These values are for extreme conditions 
and a simple optical model. For typical wing targets, the shift 
will be less than 1 y m when viewed by one tracker and less than 
2 ym when viewed by a stereo pair due to the reduced angle of view 
and smaller boundary layer thickness. 


The multiple layer model of Figure 8 was programmed. 
In order to ensure accuracy, all ray tracing variables are double 
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percision variables. The product of the index and the sine of the 
angle is computed from the length calculations and never converted 
back to an angle. This prohibits buildup of errors due to use of 
the trigometric functions. Several versions of the program were 
developed for various approximations of the index gradient. 

Figure 11 shows the index of refraction in the boun- 
dary layer for several different approximations. The zero order 
approximation assumes the index in the entire boundary layer is 
the minimum index value. Thus, a ray bends sharply at the boun- 
dary edge and is straight from then on. The first order approxi- 
mation assumes a first order polynomial, or linear variation for 
the index. A ray curve goes through the boundary layer and inter- 
sects the aircraft surface with only 40 percent of the shift indi- 
cated by the zero order trace. The second order ray trace assumes 
a second order polynomial, or parabola, for the index variation. 
This is the equation used by Shames^ to represent the velocity in 
the layer. This model leads to about a 10 percent boundary layer 
thickness error. Using this formula, the ray trace yields a 
shift of only 25 percent of the zero order ray trace. 

The actual index variation depends on the density, not 
the velocity. Using Bernoulli's equation for compressible flow, 
the density can be approximated. The Bernoulli equation is; 

y p ’ " ' ' “ ‘ ' ■ ■ ~ 

y _ 1 ^ 1“ = constant 

In boundary layer, the pressure is practically inde- 
pendent of distance from the surface and can be assumed a constant 
(p = constant). Thus: 

P = ^ 

(C^ - V ) 


^ I. H. Shames, Mechanics of Fluids , McGraw-Hill, 1962. 
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where and C 2 are constants. The boundary conditions normalized 
to free stream conditions are: 

for V = 1.0 (free stream) p =1 

V = 0 (surface) p = 0.747 

V = 1.0 = 0 

dv 

Thus: C]^ = 2.955 

C 2 = 3.955 

The final relationship between position (Y) and index 
of refraction (n) can be found by combining the Bernoulli approxi- 
mation, the parabolic velocity approximation, and the Gladstone - 
Dale equation. 

n = KP +1 
2.955 

where p = 2 

3.955 - V 


where 


V = 2Y - y2 


Figure 12 shows the path of a ray through the boundary 
layer for each model of index. The worst-case shift for a change 
from +95°C at 1 atmosphere to -195°C at 9 atmospheres can be com- 
puted. The constant index model of figure 10 yields a 7 ym shift. 
The linear index model yields a 4.2 ym shift. The parabolic index 
model yields a 3.0 ym shift. The Bernoulli index model yields a 
4.3 ym shift. These shifts are along the surface when viewed by a 
single tracker. When viewed by a stereo pair of trackers as shown 
in Figure 10, the shift along the surface is approximately one 
fifth as large but there is a shift into the surface of approxima- 
tely one-half the single tracker value. These shifts are small 
when compared to the 150 ym allowable error for the system. Thus, 
laminar boundary layers will not degrade the system. 
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Boundary layers on the tunnel window will cause an 
insignificant shift in the focus. The analysis of the tunnel 
window's effect on focus can be applied to each layer of figure 
8. The thickness and the index change in each layer is small 
resulting in an insignificant focus shift. 

Turbulent Boundary Layers 

Turbulent boundary layers may develop on the model and 
the tunnel window. The influence of turbulent boundary layers can 
be calculated by the following equation:^ 

0 = K B p /t~ 

where 

0 is the resolution limit 
K is 2.47 rad ft^/^/slug 
P is free stream air density 
t is boundary layer thickness 
B is scattering parameter 

The boundary layer thickness can be calculated by: 

0 . 37x 
^ R 

N 

where x = 1 x_ 10^ mm (typical fuselage length) 

Rjq = 120 X 10^ (tunnel maximum) 

The value of B can be found from wind tunnel 
measurements : ^ 

^ = 0.04 for M = 1.2 

pt 


® Wolfe, W. : Handbook of Military Infrared Technology. 1965. 
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Thus, the thickness (t) is 8.9 mm or 0.03 feet. The 
value of 0 is then 1.23 x 10“^ radians. This indicates that for 
collimated light, 85 percent of the total energy will be imaged 
into a spot 1.23 x 10“^ radians in diameter. In this situation, 
the light is not collimated and the boundary layer can be at the 
target or at the tunnel window. 

If the boundary layer on the tunnel window is tur- 
bulent, then the expected resolution on the model 1000 mm distant 
is approximately 1 mm. This is twice the target size and con- 
siderably larger than the 150 pm permissible position error. 
Fortunately, the system does not require that the target is 
resolved to 150 pm accuracy, but requires that the centroid of its 
image is located to that accuracy. The tracking process of 
sampling the scene brightness in a converging measurement sequence 
locates the target accurately even though the boundary layer 
precludes sharp spatial resolution. The loss in resolution will 
slow the convergence when near the final location by changing the 
slope of the error detection curve. The larger image will also 
tend to linearize the detection curve thus speeding up the first 
few cycles of the convergence. The absolute image size expected 
is quite acceptable. 

If the boundary layer is on the model, it produces a 
random target shift of approximately 10 pm. This is well below 
the 150 ym allowable error and should cause no problems. There 
may be a slight fluctuation in the intensity due to the angular 
change in the viewing of the target, but this is expected to be 
insignificant due to the extremely small angular change (0.07 
degree) and the 100 ps time scale of the tracker intensity 
measurement process. 
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Shock Wave 




Shock waves can be modeled as a single refracting sur- 
face with the index of refraction changing at the shock wave. The 
relationship between the density of the mdeium prior to and after 
a shock is:®'^ 

2 _ (y + 1) sin^ 0 
^1 24-(Y-l)M^sin 0 

where: 

= density after the shock 

P 2 = density before the shock 
M = Mach number 
0 = angle of shock to air flow 
Y = 1.47 for nitrogen at -195°C, 9 atm 

The maximum density ratio is 1.33 for a normal shock 
(0 = 90®) and M = 1.2, which is the tunnel maximum. This yields 
n 2 = 1.00937 ahead of the shock and n^^ = 1.01246 behind the shock. 
For a ray incidence of 45®, the deviation will be 0.175 degrees 
causing a target of 1000 mm distance to shift 3.05 mm. 

Fortunately, this worst-case condition will not exist for any 
reasonable aircraft model. 

A much more realistic situation is for a shock at 45® 
to the flow. This yields density ratios of 1.015 and nj^ of 
1.00923 behind the shock. For a ray incident at 45®, the 
deviation will, be 0.008 degree causing a target at 1000 mm to 
shift 138 ym. Any reduction in the angle of either the shock wave 
itself or the incidence of the light ray will reduce this value. 
The value of 1000 (range) is dependent on the actual location of 
the target. For typical wing mounted targets, a range of 700 mm 
is more reasonable, yielding a shift of only 97 ym. 

If the shock is conical., then it acts as a lens and 


^ Domasch, D. : Airplane Aerodynamics. Pitman Publishing Corp., 

1967. 
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causes a focus shift, 
be calculated by:l® 


The effective focal length of the shock can 


n„ tan 0 
2 



where : 

f = focal length 
X = distance from apex of shock 
n^^ = index after shock 
r \2 - index before shock 
0 = shock angle 


length 


For a 45° shock at M = 1.2 and x = 700 mm, the focal 
is 5 X 10^ mm. The resulting shift can be calculated by: 


Shift 



Thus, for typical wing targets (x = 700 mm) the shift 
is 97 ym in focus. This is well within the tracker's ability and 
can be ignored. 


Vibrating Objects 

The objective of the SETS is to measure deformation of 
wind tunnel models. This is done by measuring the position of 
targets placed on the model. Each target is measured by two or 
more trackers and the stereo transformation pfov'ides the correc^t 
X, Y, and Z position. 


The SETS can measure target positions at a rate of 
1700 targets per second. If 50 targets were placed on the model 
then each target would be measured 34 times every second (34 x 50 
= 1700). With more targets, there would be less measurements 


Jenkins, F. : Fundamentals of Optics. McGraw-Hill Book Company, 

1957. 
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every second and, conversely, there would be more measurements 
every second with less targets. If the model is stable, a simple 
average of the measurements would yield a target's position and 
hence the model deformation. The more measurements that are 
averaged, the less system noise and the less erroneous data points 
will affect the final position. 

The big unknown in the system is stability. Care has 
been taken in the design of the NTF tunnel to provide stable con- 
ditions and tunnel operation will undoubtedly be optimized to pro- 
vide the best conditions when data is taken. However, due to the 
extreme resolution of the SETS (150 ym of motion is maximum 
allowable), even this will not provide a non-moving model. 
Vibrations from the tunnel machinery and the model interacting 
with the air will cause the targets to move. The average position 
of the target must be measured such that the average model defor- 
mation can be found. 

The exact vibration of the target can only be esti- 
mated at this time. The expected vibration spectrum is shown in 
Figure 13 and will be a combination of sine and random components. 
Low frequency motions will be due to the sting and supporting 
structure, while high frequency motions will be from the less 
massive model itself. The highest frequencies would be expected 
at the wing tips due to vibration. Random vibration will be 
generated from the turbulent boundary layers. Thus, any measure- 
ment scheme must handle sine and random components. 

The SETS measures the target position by an interac- 
tive process that converges on the center of the target. 

Vibration during the time that a single measurement is being made 
slows the convergence but has little effect on the measured posi- 
tion due to the speed of the system. H Thus a single measurement 


Hertel, Richard J.: Study of a Stereo Electro-Optical, Tracker 

System for the Measurement of Model Aeroelastic Deformations at 
the National Transonic Facility. NASA CR-159146, 1979. 
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is the instantaneous position of the target. The SETS acts like a 
strobe in that it freezes the target position at the moment Of 
sampling. By examining these frozen position measurements, the 
average target position is obtained. 

The stroboscopic concept is ideal in examining the 
SETS data. In the typical stroboscopic application, the object is 
vibrating at a much higher frequency than the sampling light 
source. This causes the objects motion to appear greatly slowed 
and in fact stopped if the two frequencies are multiples. The 
object's true motion is not measured but can be derived from the 
measurements if taken correctly. If the object has a single 
vibration frequency, sampling at two different rates is enough to 
deduce the motion. If many vibration frequencies are present then 
many more sampling rates are required.' As long as the vibration 
is periodic and unchanging, and enough data is gathered, the 
original wave form can be reconstructed. The technique allows for 
slow sampling (compared to the signal) but requires long time 
periods (compared to the signal period) of observation. This sto- 
boscopic technique is similar to that used in sampling oscillosco- 
pes . 


Most stroboscopic investigations and all sampling 
oscilloscopes are designed to reconstruct the waveform. This is 
more information than required for the SETS where only the average 
value of the waveform is needed. Thus, many of the restrictions 
mentioned earlier on sampling at many frequencies may not apply. 

An investigation of the sampling frequencies and the resulting 
errors was performed by assuming all vibration frequencies were 
possible and sampling using various techniques would be available. 

The simplest mathematical model is to assume sine wave 
vibration and uniform sampling. A 50-Hz sampling rate was chosen 
as representative of a typical SETS data run (i.e., 34 targets 
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measured at 1700 samples/s) . Two sampling techniques were 
investigated. The first method obtained its first sample at the 
zero point in the sine wave vibration and subsequent samples at 
l/50th of a second intervals. The second method obtained its 
first sample at the positive peak in the sine wave vibration, and 
subsequent samples were at l/50th of a second intervals. In both 
cases, the signal was sampled for one second and the 50 samples 
were averaged. The results are shown in Figure 14. The maximum 
errors occur when the vibration frequency is a multiple of the 
sampling rate and, thus, the two are synchronized. The minimum 
error occurs when the vibration frequency is not related to the 
sampling rate and errors less than 3 percent of the peak vibration 
motion are obtained. 

A more complex sampling method would use multiple 
sampling rates, and the results are shown in Figure 15. Sampling 
was done at the rate of 20 per second and also at the rate of 30 
per second. The curve shows similar peak errors to the single 
rate sampling technique. Peak errors occur when the vibration 
frequency is a multiple of either sampling rate and minimum errors 
occur when the vibration frequency is not related to either 
sampling rate. Multiple rates produce multiple frequencies where 
peak errors occur but have very little effect on the average value 
of the error curve. If the vibration frequencies present are 
distributed over the frequency range, there appears to be little 
advantage to either the single or the multiple sampling rate 
methods. If, however, the vibrational frequencies are con- 
centrated at particular frequencies then a sampling scheme may be 
designed to minimize the errors. 

If Gaussian noise is added to the sine wave vibration 
of the previous models, the sampling method becomes less of a fac- 
tor in the error determination. Random noise introduces an error 
that is not removed by any sampling method and tends to mask the 
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Figure 14 Sampling at 50 Times per Second 
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errors present due to the sampling and vibration frequency reso- 
nances. In the SETS, random vibration can be expected to be a 
large component of the target motion. Thus, the error in the 
measured position will be largely dependent on obtaining as many 
samples as possible, and the sampling method will be of secondary 
importance . 

If the time between samples is allowed to vary 
slightly, the overall error will be unaffected but the minimum and 
maximum error values will be brought closer together. If, instead 
of sampling at exactly l/50th of a second intervals, the sampling 
is at intervals varying randomly between l/40th and l/60th of a 
second, the error curve peaks and valleys will be reduced. This 
is a direct consequence of the vibration frequency not being an 
exact multiple of the sampling rate. The SETS will automatically 
introduce some randomness into the sample interval due to the con- 
vergence technique used to collect data. Additional randomness 
could be introduced into the data collection if desired. 

The discussion thus far has concerned vibration fre- 
quencies above the sampling frequency. For vibration below the 
sampling frequency, a different set of problems are present. For 
these low frequency vibrations several samples per cycle are 
possible, and the portion of the signal sampled can be 
reconstructed. Figure 16 st^ows a hypothetical situation where low 
frequency vibration is present. The measured position could be 
interpreted as consisting of a 0.225Hz and a 2-Hz vibration plus 
higher frequencies. The wave can be reconstructed, but due to the 
limited time of observation (1 second data run), only part of a 
cycle is present. If the data is averaged, the low frequency 
vibration (below 1 Hz) could lead to errors. These errors can be 
reduced if the data is analyzed for these low frequencies and 
corrections made. Hopefully, these low frequency vibrations will 
be small. If they are not small, and accurate corrections for 
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them are needed, the period of data collection must be lengthened 
to use data either before or after the one second data run.^2 

Target Characteristics^^ 


Targets cannot be considered without also considering 
some aspects of target illumination and. tracking method. The 
model can move in pitch, roll, and yaw over a substantial fraction 
of the tracker system field of view. Thus, the target appearance 
can change substantially. The ideal solution is a combination of 
targets, illumination, and tracking method which yields a position 
measurement accuracy independent of changes in 


1 . 

2 . 

3. 

4. 

5. 

6 . 


7. 


Target size - magnification 
Target orientation - rotation 
Target shape - distortion 

Target illumination - magnitude, uniformity 
Target contrast - stray light 
Target focus - depth of focus 
Target viewing angle - off normal view 


^^Each wind tunnel test will consist of a series of tunnel con- 
ditrons and^ model attitudes'. "The tunnel^will 'be ■ allowed- to- sta - 
bilize after each change prior to data collection. Thus, the 
operation is a series of step and dwell operations. A typical 
SETS data run will be one second of data collection during the 
stable dwell period. 

^^Hertel, J. R. : Study of a Stereo Electro-Optical Tracker System 

for the Measurement of Model Aeroelastic Deformations at the 
National Transonic Facility. NASA CR-159146, 1979. 
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The target and its illumination must also be optically efficient 
in terms of power input to the illuminator for the signal needed 
at the camera. 

Target Size 


There are two general classes of targets; point and 
finite. Point targets are not resolved as a disc image by the 
camera and optics. Finite targets have a definable size, shape 
and even structure. 

Point targets can be defined quantitatively as having 
an image whose size is less than one-third the size of the camera 
optical spread function. The section on camera optics estimated 
the radius of the spread function to be .10-20 ym plus a 0 to 8 ym 
term for depth of focus effects. Using the dimensions of Figure 2 
and a 16 mm FL lens, the smallest demagnification is about 75 at a 
range of 1200 mm. The size of the target at the model is, there- 
fore, 2 X 750 ym/3 = 500 ym or less diameter. 

Target position for point targets is measured by 
making four samples in the vicinity of the target image. The 
position of the four samples is adjusted until all have equal 
signal amplitude. 

Numerous types of finite targets can be considered. 

In general, finite targets require more image samples than point 
targets because of the range of sizes and possible orientations. 

Two simple examples of finite targets are the white 
circle in a black surround or the black and white wedge as shown 
in Figure 17. 
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Since the model can move in range and pitch angle, 
these targets in the extreme can look like those shown in figure 
18. A scanning method which is independent of these size and 
shape changes is a raster whose dimensions are greater than the 
maximum possible target size. The scanned data is then used to 
compute the center of brightness, assuming a uniform target. This 
method is time consuming if done by computer or complicated 
electronically if done by a dedicated hardware processor. 



Figure 17. Two Examples 


of 


Finite Targets 




Figure 18. Finite Targets at Large Angles 
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A less complex or time consuming method uses a cir- 
cular scan pattern. The method seeks the black-white edge 
crossings and adjusts the center of the scan circle until two 
pairs of edge crossings are each 180° apart. This method is sen- 
sitive to target rotation. The rotation must be known or 
measured before proper adjustment of the scan center occurs. 

In general, finite targets require a tracking and 
measurement method insensitive to size and shape changes. Finite 
targets are best used where some useful information is conveyed by 
their size, shape, or pattern. 

' Active versus Passive Targets 


Targets can be either active sources of light or 
passive reflectors of light. Active point targets such as light 
emitting diodes are optically very energy efficient. All the 
energy contributes to the source intensity. Fiber optic light 
pipes are a bit less efficient due to losses along the fibers and 
at the input coupling. Nonetheless, optical fibers viewed end on 
are potentially good targets. 

Passive point targets, such as diffuse painted pat- 
terns, are optically very inefficient in this installation. Only 
a^small fraction of the. incident _energy is reflected to the’ 
cameras. Passive point targets can be illuminated by a spot light 
that moves in step with the tracker as the tracker moves from 
target to target. Such a scheme would require either premapping 
of all nominal target positions versus angle of attack, or, real 
time computation of illuminator pointing angles. It is still 
necessary to illuminate an area larger than the passive target, 
since an allowance for pointing errors must be made. 

The flux onto a target at normal incidence is given by 
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where 


A = 1 X 10“® m^, the area of a point target 
f = 2.8, the focal ratio of the lens 
SNR = 30, the signal-to-noise ratio needed in a typical 
system 

M = 125, the demagnification between target and image 
K = 2.3, the noise factor of the camera 
e = 1.602 X 10“^^ coulombs/electron 
£l £ 2 = 1 , reflection and lens transmission factors 
S = 200 y A/lumen, the camera sensitivity 
A t = 20 ys, the time to make one brightness sample 


F 


B A 
o 


4tt(2.8)^(30)^(125)^(2.3) (1.6 X (1) (1) 

(200 X 10 "^) (20 X 10~^) 


F =0.13 lumen 

The target is only 1 x 10"^ m^ in area, yet the flux is 0.13 
lumens. The illuminance must therefore be 1.3 x 10^ lumen/m^. 

This is the same as direct sunlight. For flood illumination of 
the 1 m^ field of view at the model this is at least 8 kW of 
quartz-iodide lamps in housings with very efficient reflectors. 

For a moving spot of illumination using a 10 cm^ area this is 
about 250 watt of quartz-iodide illumination with" an f/2 projector 
lens . 


A corner cube reflector used as a passive point target 
can be energy efficient. A collimated beam of light to a corner 
cube is reflected back along the same path to the source. A 
corner cube target does not exhibit the depth of focus effects 
previously mentioned because the source appears to be located at 
infinity. The two major problems with corner cubes are size and 
view angle. Also, it is not clear if millimeter sized cubes do 
not work at large angles (45°) of incidence. 

In general, active point targets such as LEDS are per- 


f erred. 
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Target Contrast 


Contrast between target and surrounding model is an 
important property. A point target tracker can, in principle, 
operate on low contrast targets. In practice, the background 
signal should be eliminated to minimize the dynamic range of the 
electronics and to remove the effects of background flux changes. 


Let Fq, Fj^, and F 2 be the average values for the 
background, sample 1, and sample 2 signals. Associated with each 
sample is a noise n^f and x\ 2 , each proportional to the square 

root of their respective signals. The track scan position error 
signal, E, is : 

(F + F ) - (F + F ) 

_ 1 o 2 o _ p “ ^9 

^ S (F + F ) + (F^ + F ) " 


F, + F + 2F 
1 2 o 


1 ‘o' ■ '‘2 

It is possible to remove the 2 Fq term by making a 
measurement of the background signal. This additional measurement 
will add 20 to 60 us to the single cycle time. 


Removal of the 2 Fq term does not completely solve the 
background problem. At null position, F]^ = F 2 f the finite noise 
associated with each of the samples becomes important. The noise, 
e, about the null position is equal to: 

2 2 

(2n + n + n ) 

^ o 1 2 

e 9P9in 

F + F„ + (2n + n. +■ n ■ ) - - 

12 o 1 2 

If Fq = Fj^, 50 percent contrast, then the position noise about 
null is at least ^ times greater than with zero background. As 
Fq becomes larger than F^, a still lower contrast, then nQ, the 
noise in the background signal becomes the dominant^^ term. 


^'^Image dissectors will respond to the 120 Hz ac flicker from 
incandescent or fluorescent lights. This can increase the 
bQ noise term considerably beyond the shot noise value. 
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If Fq <. 0.1 a condition corresponding to > 90 per- 

cent contrast, background measurement and subtraction is not 
necessary. If Fq = 0.1 Fj^ , the noise about null is only 5 percent 
greater than the zero background case; an acceptable condition. 

Tracker control loop gain depends upon the 2 Fq term in 
the error signal equation. For 0 £ F^ £ 0.1 Fj^, the change in 
loop gain is less than 1 dB, an acceptable value. 

The limits on absolute background brightness and the 
needed target intensity can be estimated for a system using LED 
targets. The required signal-to-noise ratio for the difference 
between target and background is 30:1 RMS. This signal-to-noise 
ratio allows the center of the target to be located at +0.5 ym at 
the model for zero angle at attack and center of the field of 
view. 


Let Iq and I]^ be the photocurrents entering the 
aperture of the image dissector due to the background and the 
target plus background. It is required that the signal-to-noise 
ratio of Iq^ - Iq be greater than or equal to 30. Since the image 
dissector is a shot noise limited device, the noise in each of the 
two photo currents is 




amperes 


il 



At 


amperes 
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where 


e = 1.602 X 10“^^ coulombs/electron 
k = 2.5, image dissector noise factor 
At = sample time, 20 ys in this system. 


SNR 




ek^ (I, + I 
TT 1 O 
At 


These equations assume no ac modulation of the background. If 
there is ac modulation, becomes the rms value of the modula- 
tion. As calculated here, the equation gives the greatest dc 
Iq value that is tolerable. The noises are assumed to be uncorre- 
lated and to add in quadrature. 


Figure 19 is a plot of - Iq/ the required target 
photo current versus background photo current Iq for SNR = 30 and 
t = 20 ys. An additional abscissa scale on the plot shows the 
maximum background luminance at the model for 2854 K tungsten 
illumination, f/2.8 lens and a 200 x 200 ym aperture in the image 
dissector. General office illumination is on the order of 100 
lumen m“2 steradian"^. 


The extra scale on the ordinate axis gives the target 
source intensity to cause the signal photo current. This inten- 
sity is based upon a f/2.8 lens and a red LED. Red LED's produce 
1 X 10"^ watts/steradian at 60 mA current, 25°C, and normal 
incidence. 

Undoubtedly, there will be other sources of light used 
within the tunnel. These sources could be incandescent, 
fluorescent, or laser. They may even be turned on and off for 
various reasons and will probably operate from the ac line. It is 
prudent to make the tracking system as insensitive as possible to 
these sources. 
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There are several means available to achieve this end. 
First, the targets can be made the most intense source within the 
tunnel. LED's can be pulsed to intensity levels 20-100 times 
greater than their dc, 25°C ratings for 1 percent duty cycle at 
300 Hz pulse rate. In this system, the target diode would be 
turned on only when the particular target position is being 
measured. This will require a small amount of circuitry in the 
model in addition to LED controller. 

Further reduction in system susceptibility to 
background illumination can be obained by using narrowband optical 
filters on the cameras. If a filter with a passband of 50 nm cen- 
tered on the 655 nm LED emission wavelength is used, the insertion 
loss will be about -3 dB. The out-of-band transmission is usually 
less than -60 dB. 

With this kind of optical filter, background light due 
to tungsten sources would be reduced by 13 dB, cool white 
fluorescent sources by 19 dB, but only 0 to 3 dB for HeNe laser 
sources . 


Added discrimination against HeNe laser light can be 
obtained by using a narrower bandpass and accepting a greater 
insertion loss or using the less efficient green LED's. The 
problem is the emission wavelengths of red LED's, and par- 
ticularly, the high efficiency red LED's overlap the HeNe laser 
wavelength as shown in Figure 20. 

If the narrow band optical filter does not provide 
sufficient discrimination between target and background, a phase 
coherent, modulated light system could be considered. Each of the 
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target LED's would be modulated at a frequency in the range of 


HIGH EFFICIENCY 



Figure 20. LED Relative Intensity Versus Wavelength 


5 MHz. Since the image dissector tube will pass this frequency, a 
tuned amplifier, phase sensitive detector would replace the con- 
ventional preamplifier. The system would accept optical radiation 
within the optical filter passband having an amplitude modulation 
of 5 MHz and the proper phase relationship to a reference signal. 
This is a complicated approach and is an extreme but effective 
measure. This approach avoids the 20 to 60 ys time penalty asso- 
ciated with measuring the background. The 50 percent duty cycle 
% 

of the 5 MHz modulation doubles the needed target intensity. 

Angular Distribution of the Target Light 

An almost hemispherical distribution of light from the 
target is essential. 


For a flat plate at zero angle of attack, the cameras 
view targets at a 35® angle from normal incidence. As the flat 
plate changes angle of attack, the 35° angle sweeps over a range 
of angles from 16° to 54°. The camera has a rather wide field of 
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view, 34°, therefore a nominal 35° viewing angle at the center of 

the field of view becomes 18° to 52° at the sides of the field of 

view for zero angle of attack and -1° to 71° for 19° angle of 

attack. The greatest viewing angle is about 72° for 19° angle of 

attack and occurs in the corner of the field of view. If targets 
are mounted on curved surfaces, this can lead to still larger 
viewing angles. 

For a target having a Lambertian distribution, the 
intensity varies as cos 0 , where 6 is the viewing angle from the 
normal. At e = 70°, the target intensity is only one-third of the 
normal value. Standard encapsulated LED's are not Lambertian. 
Intensity at 70° is about 10 to 20 percent of on-axis intensity. 
Stated differently, the brightness required for off axis viewing 
will determine the target source intensity. This fact alone, will 
likely dictate pulsed operation of the LED targets. 

Tracker Packaging 

The packaging concept for an individual tracker is 

f 

shown in Figure ^4^* All the optical and electrical portions will 
be maintained at room temperature and atmospheric pressure. The 
concept pictured should allow an insulated pressure vessel of ade- 
quate size to be installed in the beams surrounding the test 
-section... 


The optical system must be a custom design to fit the 
very restricted s^pace inside the beam. The concept of folding the 
optical path should provide the 40 degree field of view needed yet 
fit within the space available. 

The housing will be made from stainless steel with 50 
mm of insulation. This should provide adequate isolation from the 
tunnel pressure and temperature. The tunnel's dry air supply 
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along with electric heaters will be controlled to assure that 
internal temperatures are stable. 

Vibration of the tracker will be minimized by energy 
absorbing mounts. Analysis of the windows indicates that tilts 
must be less than half degree which appears to be easily achieved. 
Motion of the tracker affects the system in different amounts 
depending on the motion. Vibration in the plane of the window 
will simply translate the tracker and will be a very small effect. 
Vibration perpendicular to the window will change the tracker to 
object distance, but again it will be a very small effect. 
Vibration perpendicular to the window will change the tracker to 
object distance but again it will be a very small effect. 
Rotational type vibrations could cause noticeable effects espe- 
cially if the point of nutation is near the lens system. Without 
a structural design to model and a good vibration spectrum, it is 
not possible to rule out tracker vibration problems. A structural 
design and analysis is beyond the scope of this report. During 
the actual design of the SETS, a detailed mechanical model and a 
vibrational analysis must be performed to ensure that no vibra- 
tional problems exist. 

The. location of the individual trackers in the tunnel 
is shown in figures.-!', .2, and 3\ Because of the tunnel geometry, 
the_tr,acke_rs._wilJ,_n^ot _ be J.d^tical ^ _.Tr^ker^ _3 an d 4 will _vj^e^ _ 
the model at a steeper angle than trackers, 1, 2, 5, and 6. Even 
though the trackers are not identical, it is anticipated that most 
of the components can be made interchangeable. 
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RAY TRACE PROGRAM 


General Program Organization 


A series of computer programs have been written that 
trace rays in a wind tunnel. The programs are organized into 
three groups; initialization, tracing, and analysis. The ini- 
tialization and analysis programs are specifically tailored to the 
Stereo Electro-Optical Tracker System of the NTF wind tunnel and 
are limited in scope by our computer capacity. The tracing 
program is not limited and can trace virtually any ray in any 
index field with or without shock waves. 

The initialization programs load information on shock 
waves, ray source locations, tracker locations, and index of 
refraction values into various data files. While most of these 
programs are general routines and capable of handling any 
geometry, the index of refraction generating program cannot. This 
program can handle plane and conical shock waves only. Computer 
programs are available that compute the density for complex-shaped 
shock waves. These can be converted to produce the index field 
with much greater accuracy and faster than writing our own 
programs. Thus, the initialization programs were kept simple 
enough to avoid the complexities of. the large aerodynamic programs 
that already exist, yet were versatile enough to allow meaningful 
ray tracing. 


The tracing programs take the information generated by 
the initialization programs and trace rays through the wind tunnel 
model. The result of the ray trace is put into a data file for 
examination by the analysis programs. The ray tracing is per- 
formed by analyzing the medium surrounding the current ray loca- 
tion and stepping the ray to a new location. The process is then 
repeated as often as required. Subroutines handle shock waves. 
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optical elements, and all index calculations. The ray tracing 
routines will trace rays through any medium with virtually any 
shock wave shape. The restrictions described previously for the 
initialization programs are not restrictions for the ray trace 
programs. Virtually any ray, in almost any medium, with almost 
any shock wave configuration, can be traced very quickly. 

The analysis programs are designed to examine the out- 
put of the ray trace. Several analysis routines are possible, but 
only one routine is complete. That routine traces rays through 
the tracker optics and describes the image location. An addi- 
tional routine simply displays the ray trace results without the 
tracker optics. Several routines were envisioned that would 
further analyze the output and generate a new data file for a 
second pass through the ray trace routine. These two pass systems 
would provide much more accurate answers when large ray deviations 
were present. 

The general program organization is shown in Figure 
21; Software Module Diagram and Data Flow. Three main programs, 
indicated by circles, provide the basic framework. Data is pro- 
vided to and from subroutines as shown by the thin arrows. The 
initialization program and the main ray trace program build data 
files as their output as shown by the large arrows. The use of 
d^ta files allows the user to interrupt the normal ray trace flow 
to either change data values or just check them. 

Figure 22 shows the geometry used for all the computer 

programs . 

Initialization Programs 


The main ray trace program utilizes data stored in 
several different data files. These files must be loaded with the 
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Figure 21 Software Module Diagram and Data Flows 
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Figure ,'22’' Ray Trace Geometry 



f '7 L, 6 ^/ b o9‘ /<L- 


proper data for the ray trace to be meaningful. The ray ' trace 
program does not check the input data for its scientific validity. 
Thus, the operator must load all the data files correctly. The 
computer programs described in this section are designed to assist 
the operator in loading and checking these input data files. 

There are three data files for data input. Thus, there are three 
groups of programs. 

The first data file contains the ray source infor- 
mation. This file, designated as ”S" file, contains the location 
(X, Y and Z) and initial direction (direction cosines) of each 
ray. There can be from one to nine rays loaded into the file. 
Program SBLD will load this file. The program allows the user to 
input individual rays or to select one of the standard ray pat- 
terns. The standard patterns are a nine-point aircraft, shown in 
Figure 23 and a seven-point star, shown in Figure 24. The program 
asks for the initial ray location and automatically calculates the 
direction cosines and loads it into the data file. The parameters 
of the particular tracker being used are also loaded into the 
file. Any of the six trackers shown in Figure 1 can be selected 
by number. All their parameters are stored in the program. 

After the source file has been loaded, the contents 
can be examined by using program SPNT. This program prints the 

CQntenJts_ of f_ileji_ . ..The values _fpr _nine^ J?y_s_.?n^^ne _tracke^_ 

are listed. Any ray with a non-zero code number will be traced. 
The codes presently used are listed in Table 8'. The tracker 
number is listed as its code number. 

The second data file contains shock wave information. 
This file, designated as "W” file, contains the parameters 
describing all the shock waves and the optical elements. There 
are ten possible entries. The first seven are reserved for shock 
waves. Item eight and nine are automatically set to be the tunnel 
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TABLE 8 SHOCK WAVE REPRESENTATION 


Cone Code = 1 

(X - q)^ ^ (V - h)^ _ Iz - k)^ ^ Q 

C^ 

if A = B then cross section is circle 
vertex at g, h, k 

perpendicular is: T SL i; ^ T k 

A B C 

Flat Sheet (Plane) Code = 2 

A(x - g) + B(y - h) + C(z - k) = 0 
point on plane = g, h, k 
perpendicular is : Ai : Bj : Ck 


Sphere Code = 3 

(X - g)^ + (y - h)^ + (2 - k)2 - A^ = 0 
radius is A 
center of g, h, k 

perpendicular is: 2(x - g)i: 2 (y - h)j: 2(z - k)k 


Elliptic Paraboloid Code = 4 

(X - q)^ (y - h)^ (z - k) _ - 

»2 ■■ T C 

A B 

if A = B then cross section is circle 
vertex at g, h, k 


perpendicular is: Z 3 - L I: j : I K L 


Hyperboloid Code .= 5 


- g)^ _ (y - h)^ -1=0 
C^ A^ B^ 

Center at g, h, k 

distance from vertext to center = c 
If A = B then cross section is circle 


perpendicular is 


-2(x - g) -2(y - h) 2(2 - k) 

. *5 * • *> J • n" ' 


B 


Ellipsoid 


Code = 6 

(x - q)^ (y - h) 

2 2 
A B 

center at g, h, k 
Se^axes of A, B, 
perpendicular is : 




-1 = 0 




k 
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window. Item ten is reserved for the tracker entrance aperture 
and will be automatically loaded during the actual ray tracing. 

All shock waves and optical elements are described by a code 
number, three centering parameters, and three shape parameters. 

The code number identifies the general shape while the other six 
parameters give the exact mathematical formula. Any shock with a 
zero code number will not be used for the ray trace. i 

Table ^8 gives the equations and parameters for each 
shock type. This table also gives the formula for the normal vec- 
tor to the shock. This vector must be calculable whenever a ray 
intersects a shock wave. Using the geometry of Figure 22, the 
shock wave parameters must be hand calculated to best fit the 
actual shock being modeled. Once these values are known, program 
"WELD" will allow the operator to load them into the data file. 

Once the shock waves have been loaded into the data 
file, program WENT will print out the contents of the data file. 
This allows the user to check the entries and provides a record 
without changing the contents of the file. 

Program WTST tests the validity of the shock wave 
equations by plotting a cross section picture of the wind tunnel. 
Five cross sections are printed showing what portion of the medium 
has. been affected by jeach shock. This program provides a check on 
each shock wave's shape and location. 

The third data file contains index of refraction 
information. This file, designated as "N" file, contains the 
index of refraction at various points and the tunnel conditions. 
Each index value is identified by its X, Y and Z location. The 
"N" file contains data in blocks of 50 data points. The file can 
be enlarged without limit provided disk space is available. Only 
one block is in memory at a time so the file length is not 
restricted . 
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The "N" file is loaded by program START. This program 
utilizes the shock wave data file and uses supplied tunnel con- 
ditions to generate the index field. The program calculates the 
index on a rectangular grid throughout the area shown in Figure 
22. The START program is designed to generate the index data file 
for simple shock waves only. It will handle plane and conical 
shocks only. The additional shocks listed in Table 9 procude a 
complex index of refraction pattern in the tunnel that the START 
program cannot calculate. If the additional shock shapes are 
required, the index (or density) distribution must be calculated 
(perhaps by the NASA REFRC program) and put in the "N" file. 
Subsequent programs for ray tracing can handle all the types of 
shocks listed in Table ^8'. 

. i Tracing Programs 

The ray tracing programs are written to handle vir- 
tually any problem yet be simple enough to allow quick answers for 
basic situations. The ray trace programs utilize three dimen- 
sional vector calculus thus allowing for a varying index of 
refraction and shock waves at compound angles. The ray is located 
by an X, Y, Z position and its direction of travel is direction 
cosines (a, B, y ) . Once established, the ray is traced by a 
series of steps through the tunnel. The final ray location and 
direction is stored for analysis. Multiple rays are traced (nine 
in the present version) and are analyzed by themselves or in con- 
junction with other ray trace results. 

The ray tracing is done by a single main program, 
called RAY, that calls on eight subroutines for various 
situations. The main program loads into memory the shock wave 
data file and the ray source data file. It loads into memory by 
using subroutines, the nearest index data points, computes the 
index gradient and steps the ray to its next location. The user 
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is kept informed of the ray's progress throughout its journey by 
periodic status checks. Additional information is printed out 
when shocks are crossed or an unusual situation occurs. 

Ray Calculation 


Figure 25 shows the vector geometry used for the ray 
trace. The tangent at point s is the vector T(s), and the ray 
location is defined by the vector r(s). The initial ray location 
and direction are provided by initialization programs. The index 
of refraction (n) is a function of location and is provided by 
subroutines. Thus, the differential equation for light rays can 
be solved: 






= grad n 


Thus, dr/ds is a unit vector tangent to the light path. If this 
is evaluated at a point when the ray can be traced a short 
distance and the process repeated, the value of r gives the posi- 
tion of the ray (photon) at the point of interest. The value of 
dr/ds gives its direction of travel as the point. Thus, defining a 
new variable T to represent the unit tangent to the ray path makes 
the equation easier to understand without violating any mathemati- 
cal laws. - 


letting ; 


jr - dr 
ds 


reduces the ray equation to: 


dn , dT 

T -r 1- n — =; grad n 

ds ds 
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or ; 1 , ;> dn. 

^ ~ (grad n - T — ) 

ds n ds 


The term dn/ds is the directional derivative of n in the direction 
of s. The direction of s is given by dr/ds. Thus, a unit vector 
in the s direction can be calculated. 


^ = grad n A 


where A is a unit vector in the s direction, 


^ dr 
ds 


Thus, 


= grad n ' T which results in; 


^ =-i (grad' n T' (grad n • T) ) 
ds n 

Where all these terms are evaluated at point s. 


Using this result, the tangent to the ray can be found by stepping 
along the ray: 

T(s + As) = T(s) + As 


This equation can be thought of as the first terms of the Taylor 
series about a point; 

T(s + As) = T(s) + hs + h As^ 

ds ^r2 

ds ' 


or it can be thought of as the definition of a differential; 


dT _ T(s + As) - T(s) 


which reduces to the same formula. 
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The location of the light ray can be found using its 
Taylor series expansion: 

-V ?-> 

2 


2 "^ 

, A \ ■^/ A , dr(s) . , , d r(s) 

r(s + As) = r(s) + — As + ^ 


ds 


or using previously defined symbols: 


As 


ds 


-> -> dT ( s ) 

r(s + As) = r(s) + T(s) As + — - — 

ds 


As 


Thus, the ray location and direction can be calculated for each 
step . 


Data Selection 

Index of refraction data is stored on the disk in a 
single large data file. This file consists of blocks of 50 data 
points. Each data point consists of an X, Y, and Z coordinate and 
the index value. This file is read block by block into memory. 
Each block is checked for valid data points then discarded and the 
next block read in. This process provides a quick look at all the 
index data without using large amounts of memory. 

As the index data file is read and checked, the 50 
data points nearest the present ray location are retained. These 
are stored in common and are then arranged in order to distance by 
a subroutine. Thus, at any time the nearest data points are 
available. 


Gradient Calculations 


The index of refraction is a function of X, Y and Z. 
Thus, function will be approximated by a simple polynomial with 
four unknown coefficients. With four data points, the four 
unknowns can he solved for. The polynomial used is: 

N = Ax + By + Cz + D 
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This results in a linear approximation in all three directions. 

The gradient can be calculated by: 

Grad N = Ai + Bj + Ck 

A subroutine calculates the polynomial coefficients 
when given four data points. The routine converges to an answer 
only when the four data points are not in the same plane. The 
four points closest to the ray may not provide a solution. When 
this happens, data points farther away are selected according to 
Table 9. If after typing all sixteen combinations in the table 
and a solution is not reached, the index is approximated as the 
average of the four closest points and the gradient set to zero. 
This steps the ray without any direction change and provides 
another attempt at solving the polynomial. 

Boundary Representation 

All boundaries, shock waves, and optical elements are 
represented by a mathematical formula. These equations are writ- 
ten in homogeneous form where f(x, y, z) = 0. The program provi- 
des for the basic conic sections listed in Table 8. Any shape 
that can be described by six coefficients could be modeled with 
only slight modification to the subroutines. 

Each time the ray is stepped, the equations for each 
shock wave are evaluated. If it is zero, the ray is on the shock. 
If it is not zero, then its sign indicates what side of the boun- 
dary the ray is on. If the current ray step is across a shock, a 
subroutine is called and the exact crossing point is located. An 
interactive technique is used to find the crossing point. The 
equations are evaluated at a point half way between the old and 
new ray locations to see what side of the shock that point is on, 
and a new location is computed. The process is repeated until the 
crossing is found. 
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Solution 

Attempt 


First 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 


Table 9. Data Point Selection 


Data Points Used 
1 2 3 4 5 6 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


X 

X 

X 

X 

X 

X 


X 

X 


X 

X 

X 
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X 

X 
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X 

X 

X 

X 


X 


X X 

X X 


X 
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X 


X 


X 

X 

X 
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Refracting Surfaces 


Shock waves and optical elements are modeled as single 
refracting surfaces with the index of refraction changing at the 
surface. The index is calculated by the gradient subroutine, and 
the incident angle is computed geometrically from the ray direc- 
tion and the surface normal. Thus, the refracted ray angle can be 
calculated by Snell's law: n sin i = n' sin i' 

This formula holds for three dimensions but is dif- 
ficult to implement. Using the method given in Modern Geometrical 
Optics, by M. Her 2 berger, a vector solution is derived and sum- 
marized below. 


The vector form of Snell's law is: 
s X 0 = s ' X O' 

where 

s = incoming ray, direction and magnitude 
s' = outgoing ray, direction and magnitude 
0 = vector of unit length normal to shock 
X = vector cross product 

The ray vector (s, s') can be thought of as unit vectors in the 
ray direction times the index of refraction. Thus: 

s X 0 = n sin i N 

where 

N is normal to the plane s x 0. 

Thus, the original formula results in 
n sin i N = n' sin i' N 
Rewriting the original formula: 
s X 0 = s ' X 0 ' 

( s ' - s ) X 0 = 0 

^ ->■ 

s' - s = K 0 

where K is the appropriate scalar . value . 
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Solving for K; 


K = 


->• -}- 
0 s' 


- ? 


From the definition of dot product: 


s' = n' cos i' 


Substitute trigonometry identity cos x = A. - sin^ 


= n' cos 


i' = n' /l - sin^ i 


= - n'2 sin^ i' 


Substitute n' sin i' = n sin i 


= sin^ i 

Substitute trigonometry identity: sin x = /l - cos^ 

= /n'^ - n^ + n^ cos^ i 


Substitute 5 s = n cos i 


= 


n '2 n 2 + (0 1)2 


Thus, K can be calculated by: 


K = 51' - 5 


K 


= /n'2 - n2 + { 0 s ) 2 ' - 0 s 


and the outgoing ray found by; 


^ -*■ 
s' = s + K 0 
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Analysis Programs 


All analysis and printout programs read the data file 
generated by the ray trace program and print an analysis of the 
ray trace. The data file includes the ray starting and stopping 
locations, the shock wave parameters, and the general tunnel 
parameters . 


Nine rays are traced during the execution of the ray 
trace program. -This allows the analysis program to compute 
meaningful items without repeated runs of the program. In the 
initialization programs, options are built in for standard target 
locations as well as individual target selection. The standard 
targets represent an aircraft model and a star. These two pat- 
terns should allow quick analysis of most tunnel conditions. 

The analysis program LOOP traces rays through the 
tracker optics and describes the image location. It also prints 
the tunnel conditions and other useful parameters. If time had 
permitted, this program would have been expanded to generate a new 
set of ray starting conditions for a second pass through the ray 
trace routine. Thus, the name LOOP. The second pass would pro- 
vide more accurate answers when large ray deviations were present 
by forcing the rays to enter the center of the tracker's entrance 
aperture. The LOOP programs geometric optics approximations to 
calculate the ray path to the image plane. 

In geometric optics, parallel rays incident on a lens 
are brought to a focus at the focal plane. Thus, any ray hitting 
the lens is just as good as any other provided they are parallel. 
If the chief ray is selected, the ray trace is trivial. The chief 
ray is the ray that goes through the center of the lens. In 
geometric optics, the chief ray is undeviated when passing through 
the lens. Figure 26 pictures the chief ray concept. This can be 
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thought of as a pinhole camera in that the rays are undeviated. 
However, the geometric optics approximation allows a ray hitting 
the lens off center to be traced. The farther off center, the 
more error introduced by the geometric approximation. Thus, the 
computer program may be required to retrace a ray to keep its off 
centering within some limit. 

Figure 27 shows the elements needed, to derive the ray 
trace equations. The center of the lens is located at Xq , Yq, and 
Zq. The optical axis has direction numbers of A, B, and C. The 
image plane is perpendicular to the. optical axis and through the 
focal point. The focal point is D, E, and F and is given by: 

D = Xq + AK 

E = Yq + BK 

F = Zq + CK 

where K is a distance parameter found from the focal length. If 
A, B, and C hapen to be direction cosines rather than direction 
numbers, then K is the distance from the lens center to the image 
place. This distance is the focal length if the optics are 
focused at infinity. Since the restriction of making them cosines 
has not been applied, the value of K must be found. 

Focal length = f = / (Xg - 0)2 + (y^ - e)2 + (Zq - F)2 

f = / (AK) 2 + (bk) 2 + (CK) 2 

f = / (a2 + b 2 + c2) 

f = K /a2 4- b2 + 


r~2 2 2 

V A + B + C 
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The general equation for a plane through a point with coordinates 
Xp, Yp, and Zp and perpendicular to a line with direction numbers 
a, b, and c is: 

a(X - Xp) + b(Y - Yp) + c(Z - Zp) = 0 
The equation of the image plane can be written as: 

A(X - D) + B(Y - E) + C(Z - F) = 0 
or it can be expressed as: 

AX + BY + CZ = G 
where G is a constant. 

G = A(Xq + AK) + B(Yq + BK) + C(Zq + CK) 

In the geometric optics approximation, the chief ray 
is assumed to pass through the lens center traveling in the same 
direction as the ray just traced through the wind tunnel. Any 
point on this ray can be found by: 

X = Xq + ctJi 
Y = Yq + 

Z = Zq + 

where Xq, Yq, and Zq are the lens center coordinates , a , B,y ate 
V the ray direction numbers, and is a distance parameter. If a, 
3, Yhappen to be direction cosines, then % is the distance. The 
point is away from the lens center. However, the direction num- 
bers cannot be assumed to be direction cosines. 


The point where the chief ray intersects the image 
plane can be found by combining the two equations to yield: 
A(Xq + a£) + B(Yq + 3£) + C(Zq + yZ) = G 
this reduces to: 


£ 


G - AX - BY - CZ 
o o o 

Aa + B3 + Cy 


The numerator is dependent on the tracker parameters only and can 
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be calculated ahead of time. Thus: 


Aa + Be + Cy 

The value of H can be found by substituting the value of G and 
solving to yield: 

H = K(a 2 + + c^) 

thus : 

2 2 2 
^ K(A + + C^) 

Aa + Be + Cy 


It is interesting to note that the denominator is the 
dot product of the optical axis, and the chief ray and the numera- 
tor goes to unity if A, B, and C are direction cosines. 


i 
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CHARGE INJECTION DEVICES 


Real time measurement of aeroelastic deformation of 
the NTF wind tunnel models is complicated by the likely high level 
of mechanical vibration of the model and trackers. The concern 
from a measurement viewpoint is not so much with the structural 
design of the measurement system as it is with finding a measure- 
ment scheme whose accuracy and precision have low susceptibility 
to vibration. 


The system concept presented in the earlier 
report^^ used a variation on an image dissector star tracker to 
measure the image plane coordinates of self luminous targets. 
Multiple targets are tracked in seeming real time by rapid 
measurement of each target image and by achieving a capability to 
sequentially measure a large number of targets per second. This 
method requires that the stereo cameras measure the same target 
image at the same time. Key factors affecting the performance of 
the system are the amplitudes and frequencies of the vibration, 
the number of targets which can be measured per second, and the 
total number of targets. 

The vibration effects can be divided into two cate- 
gories. Category I vibrations are those which cause the ensemble 
of targets to move within the system field of view. Examples of 


Hertel, Richard J. ; Study of a Stereo Electro-Optical Tracker 
System for the Measurement of Model Aeroelastic Deformations at 
the National Transonic Facility. NASA CR-159146, 1979. 
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this vibration are the rigid body motion of the model and sting, 
the bending modes of the sting, the tracker mounting point vibra- 
tions causing image rotation or changes in the sensor line of 
sight vector. Category I vibrations can severely affect accuracy 
unless all target positions are measured in a time which is short 
compared to the vibration periods. It is category I vibrations 
which drive the system sequential target rate requirements and 
target brightness requirements. 

Example of category II vibrations are the bending 
modes of the model. These vibrations affect the sensed spacings 
among the various targets. It is precisely these inter-target 
spacings which are the key to measuring the model deformation. , It 
is assumed that the desired model deformation data is the average 
of the target positions over a period of time and not the instan- 
taneous deformation. Using this assumption as a basis for a 
signal processing algorithm means obtaining a sufficient number of 
pseudo-random target position samples over a long enough period of 
time to accommodate the bending mode frequency spectrum. 

An Alternate Approach 


Category I vibration effects can, in principle, be 
substantially reduced, if not eliminated, if all target positions 
are measured or at least recorded at the same time. This line of 
reasoning generally leads to either vidicon TV systems or fast 
framing film cameras; neither are acceptable for this application. 
MOS-charge transfer devices provide an alternative sensor 
technology . 

Developments in the field of solid state, self-scanned 
sensors have progressed to the level that such devices appear 
attractive for this application. The specific sensor device con- 
sidered here is that of the charge injection device (CID) as built 
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by General Electric. Other devices such as charge-coupled devices 
(CCD) by RCA, TI or Fairchild are also worth consideration but are 
beyond the scope of the inunediate task. 

Charge Injection Devices 


These are monolithic silicon sensors presently 
available as 128 x 128 pixels of 45 x 45 ym or 244 x 248 pixels of 
36 X 46 urn. The later device is available in a camera with RS170 
compatible video and sync. A 256 x 256 pixel sensor with 20 x 
20 ym pixels is reported to be in development and a camera is 
expected in 1982. A key property of the CID is that the pixels 
are contiguous. 

Overall System Description ^Using CID 


The CID sensors would use the same tunnel locations 
and generally the same optical concepts as described in this and 
the previous report. By a combination of optical bandpass 
filters, target brightness, and lens f/number, the, targets would 
appear as bright spots in an otherwise dark field. The optics 
will be adjusted deliberately off focus so that the otherwise 
small image of the target covers 9 pixels on the CID. This multi- 
pixel coverage is needed to implement the tracking algorithm. 
Multi-pixel coverage restricts the minimum inter-target distance 
for all model orientations and camera angles. If two target ima- 
ges merge, the signal processing algorithm becomes much more 
complex and may not yield a solution since targets are viewed 
simultaneously. 

A coarse degree of real time tracking is needed in 
order to minimize the volume and rate of data sent to the computer. 
The CID tracker extracts the nine pixels associated with each 
target image from the full frame and only transmits these nine. 
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Table 10 lists the system performance objectives as 
they relate to a CID based system. 

System Performance Analysis 


The performance will be examined for three different 
density CID's: 

• 128 X 128 pixels of 45 x 45 ym 

Available now 

• 256 X 256 pixels of 20 x 20 ym 

Available early 1982 

• 512 X 512 pixels of 20 x 20 ym 

Projected 1984. 


CID performance data will be that of the 128 x 128 
device scaled to that of the other densities. The results are 
summarized in Table 12 ''. (/ 

Field of View, Lens Focal Lengthy and Magnification 


The nominal distance from model to tracker is about 
1600 mm. For simplicity, a 1 x 1 m field of view will be assumed 
to exist at this distance. 

The required lens focal length is 

. f £ _ image size 

1600 field of view 


= Magnification 




Table 10. SETS Characteristics - CID Version 


Sensed Volume 
Target Sources, Number 
, Type 
, Size 

Model Deformation at 
Target Location 

Sensor Position Deviation 
Detection 

Stereo Position Deviation 
Calculation 

Model Vibration Limits 


0.8 m by 1.4m by 0.5 m 
5 to 50 

Light-Emitting Diode 
0.5 mm diameter 


+ 20 mm maximum 


0.05 mm all axes 


0.1 mm all axes 

+ 12.5 mm at 20 Hz 

+ 5 mm at 50 Hz 

+ 1.25 mm at 200 Hz 


Sensor Environment Control 

Tracking Electronics 

Calibration, geometric 
Pre-Test Reference 
Data Processor Type 
Data Storage 
Control/Display 
Position Display 
Hard Copy 


70°F Dry Air at 130 psi (NASA 
Supplied) Electrical Heating 

Digital, in Tracker head 
Software, in computer 

Area map, prior to test 

Targets on Stable Portion of model 

MODCOMP 

Magnetic Tape 

Keyboard/Printer 

CRT Display 

Line Printer/Plotter 
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Sensor Image 

( Pixel ) Size 


128 

X 128 

5. 

76 

X 

(45 

ym) 

5. 

76 

mm 

256 

X 256 

5. 

12 

X 

( 20 

ym) 

5. 

12 

mm 

512 

X 512 

10 

.2 

X 

( 20 

ym) 

10 

.2 

mm 


Focal 


Length 

1/M 

9.2 mm 

174 

8.2 mm 

195 

16 . 4 mm 

98 


Target Position Determination, 


The objective performance is to resolve 50 wm change 
in position at the model. This dimension can be referred to the 
tracker sensor as (50 ym)(M) = d. The value of d for the respec- 
tive cameras are 0.3, 0.25 and 0.50 ym. 


The plan is to use a centroid finding algorithm. ^2 
This method uses the numerical intensity data from the 3x3 array 
of, illuminated pixels to calculate the following: 


11 

12 

13 

21 

22 

23 

- 31- 

32- 

- 33 - 


Salomon, P.M. : Charge Coupled Device (CCD) Tracker for High- 

Accuracy Guidance Applications. Optical Engineering, Jan/Feb., 
1981, Vol. 20 No. 1, p. 135. 
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s = i: E 

i=l j=l 


The centroid offset from cell 22 is found by 



Rf - 

y = ky S 

when Kx and ky are scale factors. 

- - - This algorithm, in effect, - must -inf er the center, of 

the target image to a small fraction of the pixel size. 


Sensor 

Pixel 

d 

Pixel/d 

128 

45 y m 

0.3 

150 

256 

20 M m 

0.25 

50 

512 

20 y m 

0.50 

40 
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The column pixel/d is a measure of the signal to noise 
ratio required to accomplish the desired measurement. Data sheets 
for the present 128 x 128 CID give spatial SNR = 25 typical and 
temporal SNR = 256 typical. These are presumably, for full scale 
signals and are likely proportional to the signal. Observation of 
two 128 X 128 CID cameras would say that the 256 value is 
optimistic; 64 is closer to an observed, real world value. 

The spatial signal-to-noise ratio is the pixel-to- 
pixel variation in sensitivity and dark current. These variations 
are presumed to be linearly correctable given two data bases; the 
dark current at each pixel and the normalized sensitivity at each 
pixel . 


One would conclude that the 256 x 256 element CID 
might be satisfactory but the 128 x 128 element device will be 
severely taxed if not incapable of the desired performance. The 
512 X 512 device is not expected for another two years. 

Frame Rate 


CID's are MOS semiconductors. Such devices have maxi- 
mum clock rates of the order of 1 x 10^ Hz. The 128 x 128 camera 
is rated at 500 K pixels per second data rate and uses a clock at 
5 X 10^ Hz. It is assumed that the 500 K pixel limit is a speci- 
fic choice by the manufacturer to permit 30 frame/second operation 
and not the absolute maximum rate. A second 244 x 248 TV com- 
patible device is available which is specified at 4.5 M 
pixels/seconds . 

Thus, one might expect that the maximum frame rates 
could be 275, 65, 18 frames/ second for the 128, 256, and 512 
devices . 
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The 65/second frame rate in conjunction with a pulsed 
target array will be capable of freezing the vibration motion and 
following the 65/2 Hz and lower frequency motion. 

Data Rate to the Computer 

The data rate per tracker is estimated on the folowing 

basis ; 

• 32 targets/f ield of view 

• 9 pixels/target 

• 275, 65 or 18 frames/second 

• data compression to extract only the nine (9) 
pixels/ target 

• Data transmitted 

10 bits brightness/pixel 

14, 16, 18 bit brightest pixel address/target 

The rate decreases because of the slower frame rates 
for the larger pixel count devices. If all devices are assumed to 
be capable of at least 30 f /second, a likely development objective 
for a 512 device, then the 12K rate becomes 20K. 
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RATE CALCULATION 





128 

256 

512 

Targets 


32 

32 

32 

pixels/target 

X 

9 

9 

9 

Video bits/pixel* 

X 

12 

11 

10 



3456 

3168 

2880 

targets 


32 

32 

32 

Address bitss/target 

X 

14 

16 

18 



448 

512 

576 


Bits/Frame 

3904 

3680 

3456 

Frame Rate (max) x 

275 

65 

18 

Convert to bytes v 

8 

8 

8 

by tes/sec 

134K 

30K 

8K 

Reserve @ 30% 

40 

9 

2.4 

Data Link Overhead @ 20% 

26 

6 

1.6 

Max Data Rate 

200K 

45K 

12K 

* bits/pixel = 1.44 In (6N/F) 

where F is the fraction of 

full scale 

signal 

the 


target generates. For trackers 0.25 < F < 0.7 
N is pixel/d 
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The 512 X 512 devices can likely acconunodate more than 
32 targets while the 128 x 128 may. not be able to handle even 32 
without some adjacency effects. 

The image dissector tracker, operating at an average 
1500 targets/second and a peak rate of 5000 targets/second, genera- 
tes 9 K bytes to 30 K bytes/second. 

Target Address 32 bits 

Reserve and data ,, ... 

16 bits 

Link overhead @ 50% 

48 bits = 6 bytes 

Hardware Differences Between CID and IDT Approaches 


Both sensor approaches can use the same concept of 
packaging; a temperature controlled, pressure isolation vessel. 

The circuitry, which is internal to the vessel, is quite different 
for the two cases. 

The IDT is an active tracker whose task is to maintain 
lock on multiple targets, measure their position, and transmit 
position data to the computer via a serial data link. There is 
analog and digital circuitry present, A/D and D/A converters, and 
a signal processor for the tracking function. There will be 
parallel-to-serial data buffers for the coordinate data. 

The CID approach will have mainly digital circuitry, 
and probably one or two full frame memories, to buffer the CID from 
the signal processor. The signal processing module will sense the 
brightest pixel of a group of nine and format the data for 
transmission . 
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At the computer end of the data link, so long as the 
CID and IDT data rates are compatible, there will be little hard- 
ware impact. 


Software Differences Between CID and IDT Approach 

There will be significant software differences because 
the IDT provides target coordinates, whereas the CID concept does 
not. The CID data must be processed by the interpolation 
algorithm to extract target centroid data. This process is pre- 
sently assigned to the computer and is intended to be a post pro- 
cessing task because of the number of unknowns associated with 
dark signal and pixel-to-pixel sensitivity variation corrections 
required. 


Data Link 


The hardware section cited the use of a serial data 
link for both IDT and CID. The concept is to use two separate 
lines for the data link; one is data from the computer, the other 
is data to the computer. Each tracker would have its own pair of 
lines for a total of 12 plus 2 for the LED controller. (Power 
lines are additional.) 

The data rate estimates for the IDT and 256 CID are 
both of the order of 50K bytes/sec. Harris Semi-conductor makes 
an Ayschronous Serial Manchester Code Adapter (HD-6408) capable of 
encoding and decoding high speed (lOOK byte per second) asynchro- 
nous serial data streams. A computer/tracker data link based on 
this device is illustrated in Figure 28. The clock rate 
illustrated is the maximum for the device and corresponds to a 
lOOK byte/second rate. Aside from the data source, the signals 
controlling the latch, and the start of transmission, the scheme 
illustrated will suffice for either the IDT or CID sensor 
approach. 
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Data Source 




12M Hz 


Clock 


T 

100 M 


12M Hz 
Clock 


Figure 28 

100 K byte/second Data 
Link for Either IDT 
or CID Tracker 











/ 2 O & ^ 


In either case, the computer is expected to receive 
its input via a direct memory processor. Data from the computer 
will likely be via a register I/O port and at a lower rate. The 
same link can be used for the data to the trackers or LED 
controller but likely at a slower clock rate. 

Data given in a Texas Instruments application note 
shows that a link such as this composed of SN75116 drivers and 
receivers and using Beldon 8795 wire will accommodate 2.5 M byte 
rates at 100 meters. 
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DISTORTION CORRECTION 

The SETS camera, optics and associated windows will 
introduce distortion into the SETS data. A correction for this 
distortion will be made using the raw data directly from the SETS 
before any stereo transformations are made. A correction prior to 
the stereo transformation is required because the distortions 
introduced are angular errors rather than positional errors. 

Since the SETS tracker is, in essence,. an angle measuring device, 
the correction should be applied directly to its data. 

The distortion from the optics and the camera cause 
the image to appear slightly displaced from its expected position 
on the image plane. These shifts would directly affect the 
accuracy of the system if they were not corrected. Some distor- 
tion is present in all optical and camera systems and does not 
necessarily indicate a design or construction problem. The SETS 
can easily incorporate a distortion correction method since the 
data is already digitized and stored in the computer. 

Camera and Camera Lens Effects Upon System Precision 

The camera and lens act in several ways to limit the 
system precision. The simplest limitation to understand is the 
one due to small but finite sized deflection increments. The 
demonstration system, described in the Final Report for the SETS, 
NASA Report Number 159146, used 13-bit deflection D/A converters 
over a 4.1 x 4.1 mm field of view. The minimum step size was 
0.5 m at the camera faceplate. The system design concept pre- 
sented in this report uses a 10 x 10 mm working field of view. 
Fourteen and 15-bit converters provide 0.8 to 0.4 urn increments, 
respectively, while covering a 12.5 x 12.5 mm scanable field of 
view. Such converters are available for the camera. 
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D/A converters and cameras have temperature dependen- 
cies which affect the size of the area scanned and the stability 
of the scan center within the field of view. A change in scan 
size is equivalent to a change in scale factor at the camera 
faceplate. The magnitude of the error is equal to the distance 
from the center of the field of view times a + 20 ppm/°C tem- 
perature coefficient. Values range from zero at the center to + 
0.12 Um/®C at the edge of the field of view. 

The centering term is + 10 ppm/°C of the full scale 
size of + 10 X 10“6 12.5 mm = + 0.12 ;im/°C. All coordinates 

within the field of view are equally affected by this factor. The 
pressure isolation container used to house the cameras must have 
temperature control because of the -195°C to 95°C tunnel environ- 
ment. The temperature coefficients require the control of the 
temperature, at certain key locations in the container, be the 
order of plus or minus 1°C to 2®C. 

The mechanical packaging of the camera and its lens 
must prevent shifts of the tube relative to the focal plane or 
coil system. The need is for less than + 0.5 ym instabilities, a 
value consistent with the needs of precision star trackers used 
for space flight. 

\ 

The camera and lens both introduce distortions of the 
field of view. Figure 29 shows the magnified distortion of an 
uncorrected image dissector camera without lens. It is assumed 
the lens distortion is comparable or less in magnitude. The 
dotted boxes are on a 1 mm grid and centered on the true position. 
The distance from the center to the side of the smaller boxes is 
+ 25 ym; + 40 ym for the larger boxes. 
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This distortion is, in principle, correctable using a 
lookup table and interpolating between table entries. A key 
question is, how large must the table be in order to correct 
+ 0.5 ym? 


Consider a group of four boxes in Figure 30 whose cen- 
ters are denoted by A, B, C, and D. The distance between AB is 
2 ram, and the same holds for CD. Assume the x and y position 
errors are known by measurement for the four positions. The posi- 
tion error for any point E inside the rectangle ACBD as shown in 
Figure 30 can then be calculated by interpolation using 

Ex = k^Cjj + (1 - kj^) Djj 

Ey = k2Ay + (1 - k2> By 

where the subscript denotes the component of the error and k^^ , 
k 2 is the position of point E from C and A, respectively. 

Analysis of the 157 points in Figure 29 showed E^ and 
Ey calculated by interpolation different from the measured values 
for Ej^ and Ey by + 2 ym typically and + 4 ym at the greatest. 

Assume the interpolation error is proportional to the 
area of the ACBD cell. Then a + 0.5 ym interpolation error can be 
achieved using ( 2 )^ x 0.5/4 = 0.25 mm^, 0.5 x 0.5 mm cell. The 
12.5 X 12.5 mm will be mapped to a sufficient density using 
(12.5)2/(0.25) = 625 points. 


System Calibration 


An array of approximately 625 points will be placed in 
object space where the aircraft model would be. The location of 
each point will be measured, and the correction required will be 
calculated. These measurements will be the data base for all 
distortion correction. Since the data is taken using the entire 
system, it will remove the distortions from the windows, optics 
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and camera in one step. These corrections are valid for the 
tracker until it is physically modified or adjusted. Periodic 
calibration may be required due to aging and the abusive environ- 
ment of the tunnel. 

Due to the extreme temperature and pressure range the 
tracker experiences, it may be necessary to have calibrations 
under various conditions. If this does become necessary, a data 
set can be taken for various conditions and the correct set used 
for each data run. There would be _no additional post processing 
time required, but the initial measurement of the calibration data 
sets would take more time. 

Post Processing 


Since the calibration data consists of measurements of 
specific points, some method of applying it to a, random point is 
required. Previous sections showed that linear interpolation is 
adequate and the method presented there is the simplest form of 
interpolation and is quickly computed. Even though four data 
points are used (A, B, C, and D) , only two (A and B) are used for 
the X direction and two (C and D) for the Y direction. This uti- 
lizes only half of the data available, as the y errors at A and B 
and the x errors at C and D are not used. It also makes the X and 
- Y corrections independent as. separa_te points are used for each. 

The interpolation is least accurate for the x correction when the 
point is near A or B and is least accurate for the Y correction 
when near C or D. This can be demonstrated by Figure 31. When 
point E is near point A, the X axis correction can be found using 
either CD or FG. Both pairs will predict the value at E yet only 
one pair is used. The Y value at E is known exactly as it is 
measured at point A. Thus, while this method may be the easiest 
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FIGURE 3l Point E Close to Point A 
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to compute, it does not efficiently utilize the data available and 
introduces some unwanted effects. 

A preferred method utilizes all four data points for 
both the X and y corrections. This can be done by fitting a 
polynomial to the four data points. Thus, the error in the x 
direction at a point (x, y) is; Error = + C2Y + C3XY + €4^ 

where the coefficients are evaluated for the direction. A similar 
equation provides the error in the Y direction. Using four data 
points, the four coefficients for each direction can be found. 

The matrix form of the four equations is; 


Exi 


Cl 


Xi 1 

^X2 

= 

C2 


X2 Y2 X2Y2 1 

^X 3 


C3 


X3 Y3 X3Y3 1 

^X 4 


C4 


X4 Y4 X4Y4 1 


Using matrix manipulations, the coefficients can be found and the 
polynomial used to find the error. 

The polynomial contains a linear term in X, a linear 
term in Y, and a cross term of XY. This combination generates a 
true linear approximation. The cross term does not introduce a 
higher order correction as might be expected. 

For true linear interpolation, the value is that pre- 
dicted by drawing a straight line between the data points. Figure 
32 shows the process involved. The four corner points (A, B, C, 
and D) are measured data points. From these four points, the 
error is found at four points on the edge of the square with the 
same X or Y coordinate as the final point. The error at the final 
point is then estimated both vertically and horizontally. These 
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two estimates are always the same even though they are arrived at 
by dif f erenf^means . This method can be expanded as shown in 
Figure 33^for any four random points. All of the interpolation is ^ 
the same, only the mathematics is more complex. 

The polynomial approximation provides exactdy the same 
answer as the linear interpolation method. Figure J,4'^is a graphi- 
cal representation of the linear interpolation method. The X-Y 
plane is the location of the point, and the vertical is the magni- 
tude of the error. The error is known at the four corners (A, B, 

C, and D) and is found everywhere else by using straight lines 
that are parallel to either the X or Y axis. The polynomial is 
the equation of the surface of all such approximations. This can 
be seen by holding either X or Y at a constant value: 

Error = Cj^X + C 5 + CgX + C 4 
which reduces to 

Error = + Cg)X + C 4 + C 5 

which is a straight line of the form E = mx + b. Thus, the two 
methods produce identical results, but the polynomial method has 
distinct computational advantages. 

Using the polynomial approach allows the data points 
to be randomly located without increasing the computation time. 

This allows an automated system to select points or dicard bad 
data points. If there is an area of particular concern, either 
for better data or a difficult distortion area, additional data 
points can be collected in that area. A computer program was 
written in FORTRAN IV using matrix manipulations to solve the four 
point polynomial. Each solution requires approximately 160 milli- 
seconds or about six per second. With the SETS taking data at a 
rate of 1800 targets per second, a one-second data run would 
require about 10 minutes of post processing. This time can be 
reduced drastically if the polynomial for each area is computed 
once and then stored for use should a point again fall in that 
area. 
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FIGURE 34 Four Point Linear Interpolation 



CONCLUSIONS AND RECOMMENDATIONS 


None of the error sources examined would prevent 
meaningful measurements of model deformations through the use of 
the Stereo Electro-Optical System, provided care is taken in its 
design and adequate post-processing of the data is performed. 
Detailed analysis of the actual design should minimize most of the 
errors to acceptable levels. Vibration, of the tracker may prove 
to be the most difficult area to minimize due to the complex 
nature of the problem and the uncertainty of the vibrational 
environment. Shock waves can produce large errors that are 
correctable, provided the exact location of the shock is known. 
This post processing to remove shock-introduced errors is essen- 
tial and depends on knowing the shock wave locations and shapes. 
Without this correction, or one based on assumed locations, the 
system performance will be limited. 

Table 12 lists the possible errors when viewing a 
model using the Stereo Electro-Optical System. All the areas of 
concern can either be corrected or are small and do not degrade 
system performance. All of the corrections require accurate 
information about the tunnel conditions at the time the data was 
taken. Much of this information is readily available such as tun- 
nel pressure, temperature, Mach number and Reynolds number. Some 
of the information may be available only after the tunnel data run 
has been analyzed. 

The windows on the tunnel wall and in the tracker are 
distorted due to pressure and temperature. The flat windows 
become very weak lenses. The effect is to slightly defocus the 
image. The defocus is such a slight amount that it is insignifi- 
cant. Tunnel window tilts relative to the optical axis produce a 
shift in the target location that could equal the entire error 
budget of 150 micrometers if the tilt were severe (half degree). 



TABLE 12. SUMMARY OF POSSIBLE TUNNEL INACCURACIES 
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Care in the design should prohibit severe tilting to limit the 
error introduced to well below allowable limits. 

Changes in the pressure and -temperature of the tunnel 
atmosphere have a direct impact on the target location as seen by 
the tracker. Both increasing the pressure and decreasing the tem- 
perature cause the atmosphere to become more dense and thus 
increases the optical index. This. increase causes the rays to 
bend at the windows in a non-linear fashion. The model appears to 
grow (4 percent in the worst case) but not uniformly. Some por- 
tions will grow as much as 9 percent while other sections grow 
only 1 percent. The data can be corrected to remove this effect, 
provided the tunnel conditions are known. 

Laminar boundary layers on surfaces produce a gradual 
change in the optical index. This gradient causes the light path 
to bend depending on the properties of the layer and the incident 
shift in the focus. A layer on the model will cause a target to 
appear displaced horizontally and into the model. The apparent 
displacements are much smaller than the system resolution and will 

i- 

not be a problem. 

Turbulent boundary layers on surfaces produce a fluc- 
tuating optical index causing the ray path to change randomly. If 
the layer is on the tunnel window, the fluctuations cause the 
image to blur reducing the optical resolution. The degraded reso- 
lution is still well within the design limits. If the layer is on 
the model, it produces a small random shift that will not affect 
system performance. 

Shock waves produce a sharp change in the optical 
index and act as single refracting surfaces which bend the optical 
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rays. The amount of bending and the effect on the optical ray 
depends on the shock strength, the angle of incidence between the 
shock and the ray, and the distance from the model at which the 
ray intersects the shock. Optical rays . intersecting shocks near 
the model are only slightly affected and correction is not 
necessary. Rays which intersect shocks near the tunnel window 
will be significantly affected and must be corrected. Correction 
can be made however, only if the strength, shape and location of 
the shock on the model are known. Curved, shocks produce a lens 
effect which cause the point at which the target image is focused 
to change. This focus shift may or may not be significant 
depending on the individual shock. 

Motion of the tracker or the target during a sampling 
period introduces a slight error into the convergence routine. 
Since the target is sampled several times to get its location, any 
motion during the sequence leads to errors and may require addi- 
tional samples. Based on the vibration levels predicted by NASA, 
acceptable accuracy and sample times are expected. 

Motion of the tracker or the target between sample 
periods can lead to large errors if all trackers do not sample 
simultaneously. Additional problems could develop if both the 
sampling and the vibration are periodic. They could combine to 
produce stroboscopic effects that could give errors. Low fre- 
quency vibration can also produce errors if the data collection 
times are short compared to the period of the vibration. The 
sampling technique proposed reduces all these errors to acceptable 
levels by adjusting the data collection time, varying the target 
order, and synchronizing all the trackers. 

Low target intensities can lead to problems as the 
convergence time increases. The target and background levels 
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anticipated should provide adequate signal-to-noise ratios for 
normal operation. 

The influences on tracking accuracy of the factors 
cited above appear to be definable to the extent of describing the 
general level of the effect on system output. 

The use of a 128 x 128 pixel CID-type system does not 
appear practical. If the technology advances as expected, it will 
be a realistic alternative in the future. It is possible to 
design a data link that is device independent to allow for future 
system improvements. 

The computer program developed does trace virtually 
any ray through the complex tunnel conditions. The program can 
handle shock wave, boundary layers, pressures, temperatures, 
nonhomogeneous index of refraction, and optics of any conic sec- 
tion or polynomial of sixth order or less. The program does 
approach the goal of a universal ray trace program capable of 
handling any tunnel condition or ray geometry. 

With the cautions noted to limit window tilt and other 
causes, it now appears that the typical range of effects antici- 
pated for the NTF do not preclude successful operation of the 
tracking system. 
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